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ABSTRACT 

We present 18 GHz Australia Telescope Compact Array imaging of the Mpc-scale quasar jet 
PKS 0637-752 with angular resolution ^ 0'.'58. We draw attention to a spectacular train of quasi- 
periodic knots along the inner 11" of the jet, with average separation ~ 1.1 arcsec (7.6 kpc projected). 
We consider two classes of model to explain the periodic knots: those that involve a static pattern 
through which the jet plasma travels (e.g. stationary shocks); and those that involve modulation 
of the jet engine. Interpreting the knots as re-confinement shocks implies the jet kinetic power 
Qjet ^ 10"*^ ergs s~^, but the constant knot separation along the jet is not expected in a realistic 
external density profile. For models involving modulation of the jet engine, we find that the required 
modulation period is 2 x 10'^ yr < r < 3 x 10^ yr. The lower end of this range is applicable if the jet 
remains highly relativistic on kpc-scales, as implied by the IC/CMB model of jet X-ray emission. We 
suggest that the periodic jet structure in PKS 0637-752 may be analogous to the quasi-periodic jet 
modulation seen in the microquasar GRS 1915-1-105, believed to result from limit cycle behaviour in 
an unstable accretion disk. If variations in the accretion rate are driven by a binary black hole, the 
predicted orbital radius is 0.7 ^ a ^ 30 pc, which corresponds to a maximum angular separation of 
~0.1-5 mas. 

Subject headings: galaxies: active — galaxies: jets — quasars: individual (PKS 0637-752) 



1. INTRODUCTION 

The jets of FRII radio galaxies and quasars often ex- 
hibit regions of enhanced brightness, commonly referred 
to as "knots". The production of knots and the asso- 
ciated particl e acceleration mechanism are poorly un- 
derst ood (e. g. lStawarz et aLll2004HNiemiec fc Ostrowskil 
|2006[) . Jet knots have traditionally been identified with 
strong shocks in a uniform and continuous flow. In this 
class of models, particles are accelerated at the shock via 
the first order Fermi mechanism, resulting in a power- 
law electron energy distribution. Such internal shocks 
could be formed due to intrinsic velocity irregularities 
in the jet (|Reesll 19781) . they could be re- confinement /re- 
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collimation shocks (jKomissarov fc Fallel[l998h . or they 
could be a ssociated with large scale in stabilities in the 
flow (e.g. iBicknell fc BegelmanI [T996() . Alternatively, 
knots in large-scale, powerful radio sources may repre- 
sent moving and separate portions of the jet matter with 
higher speed and/or luminosity (Stawarz ct al. 2004). 

The radio morphology of PKS 0637-752 (z = 0.653) 
consists of a bright core; a knotty, one-sided radio jet 
extending 15 arcseconds west of the core; and a counter- 
lobe 11 arcseconds east of the core (see Figure [T]). The 
apparent superluminal motion of pc-scale jet components 
wit h v/c » 13 implies a j et viewing ang le less than 
9° (|Lovell et all [2000 : Ed wards et al.ll2006[) . and hence 
source size greater than 1 Mpc. The object has received 
considerable attention in the past decade, largely due to 
the bright X-ray emission associated with the 100-kpc- 
scale western jet. The strong jet X-ray emission is hard to 
explain in terms of standard emission mechanisms such 
as t hermal Bremsstrahlun g or synchrotron self Comp- 
ton ([Schwartz et al.l I2000D . A s a po tenti al solution to 
this p roblem, iTavecchio et al.l ()2000D and iCelotti et al.l 
()2001t ) proposed the now-popular beamed, equipartition 
inverse Compton model, in which the cosmic microwave 
background (CMB) provides the soft photons that are 
scattered to X-ray energies (IC / CMB model) . 

Since the launch of Chandra in 1999, several tens of 
quasar X-ray jets similar to PKS 0637-752 have been 
detected. Despite a decade of intense observational and 
theoretical analysis, the X-ray emission mechanism, dy- 
namics, and the physical processes shaping jet mor- 
phology remain unresolved. A number of problems 
with the IC/CMB model have been identi fied (see e. g. 
[Harris fc Krawczvnskil 120061: iWorralll [20091) . and whilst 
none of the problems is seen as directly refuting the 
model, mounting criticism has meant that competing 
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Fig. 1. — : ATCA image of PKS 0637-752 at 17.7 GHz. The colour scale is logarithmic between a minimum and maximum surface 
brightness of 0.5 and 38 mjy/beam respectively. This colour scale was chosen to show clearly the regular nature of the inner jet knots. 
The restoring beam is circular with FWHM = 0'.'575. The dynamic range is 16,000:1; the off-source rms noise level is 0.3 mJy/beam, and 
the peak of the faintest knot in the inner jet is 5 mJy/beam. Knots in the inner jet are real and are well above the side-lobes from the 
bright core which are at or below the noise level. 

models such as d ual population synchrotron models (e. g. 
iTavecchio et al.|[200 3: Jester et al. 2 0061 ) are considered 
as likely candidates. Additional constraints are required 
to determine the X-ray emission mechanism and the na- 
ture of large-scale quasar jets. In the present work, we 
consider a new constraint on the jet dynamics: the exis- 
tence of quasi-periodic peaks in the jet emission. In |J2] 
we describe the radio observations and data reduction. 
In f}3] we discuss the periodic structure in the jet and 
possible interpretations. 

We adopt the following cosmology: Hq — 
71 km/s/Mpc, nm = 0.27, Qa = 0.73. 



gle to the line of sight must be greater than ^ 4° , other- 
wise the deprojected source siz e would be uncomfo rtably 
large (> 2 Mpc). Based on the iLovell et all (|2000( ) com- 
ponent identificati ons and inferred super luminal motion 
in the pc-scale jet, lEdwards et al.l ()2006[ ) obtain a limit 
oi 9 <9° for the jet viewing angle. The par sec-scale and 
arcsecond-scale radio jets a re well aligned (jLovell et al.l 
120001 : iSchwartz et al.l l2"000( ) so that, unless the jet goes 
through a large bend in a plane perpendicular to the 
plane of sky, the jet viewing angle on arcsecond scales 
is likewise < 9°. Thus, if the knots are stationary, the 
de-projected knot separation is 



2. OBSERVATIONS 

PKS 0637-752 was observed with the ATCA in the 6C 
configuration at 17.7 GHz and 20.2 GHz on 10 May 2004. 
A full 12 hour synthesis was obtained, recording 128 MHz 
bandwidth. Regular scans on the nearby phase calibrator 
0454-810 were scheduled throughout the observations, as 
well as scans on the ATCA flux calibrator PKS 1934- 
638. Standard calibration and editing procedures were 
carried out using the MIRIAD data analysis package. 
Following the initial calibration, data were exported to 
DIFMAP and several imaging/self-calibration iterations 
were performed, incorporating both phase and amplitude 
self-calibration. 

3. PERIODIC JET KNOTS 

Trains of bright knots, in some cases with a regular or 
quasi-periodic pattern, occur in radio galaxy and quasar 
jets: 3C 175, 3 C 204, 3C 263, 3C 334 and 3C3 03 are 
prime examples (jBridle et al.lll994] : lKronbergll 986^. Fig- 
ures [T] and [2] illustrate the periodic nature of the jet knots 
in PKS 0637-752, extending along the inner jet up to 11 
arcseconds from the core. The average apparent knot 
separation is I'.'l with standard deviation 0'.'15, corre- 
sponding to projected separation Dapp = 7.6 ± 1.0 kpc. 
Following .Mehta et al.. (,2009. ) . we argue that the jet an- 



50 kpc < D < 110 kpc 



(1) 

As discussed in §3.2| if the knots are discrete jet com- 
ponents with relativistic velocity then light travel time 
effects must be taken into account when calculating the 
knot separation. 

The periodic nature of the knots in PKS 0637-752 is 
striking, and here we ask the question: what physical 
process is responsible for the periodic knot structure? 
We can consider the knots to be either: (A) a static 
pattern through which the flow travels (e. g. stationary 
shocks); or (B) a real variation within the flow caused 
by a variable jet engine. Section 13.11 deals with models 
that fall under category A, while Section [3?2] deals with 
models that fall under category B. 

We defer a more detailed investigation to an upcoming 
publication using spectral index and polarisation data 
based on more sensitive, and higher resolution, ongoing 
ATCA imaging. 

3.1. Large-scale shocks in a continuous flow 
3.1.1. Re- confinement Sh,och:s 

Hydrostatic jet conflnement is accompanied b y the for- 
mation of regularly spaced shocks along the jet () Sanders! 
il983i : iKomissarov &: Falls, 1998.) , which could explain the 
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Fig. 2. — : Plot of jet brightness integrated along the North- 
South direction, as a function of distance from the core in a westerly 
direction. The bar in the top left of the plot indicates the off-source 
rms of the projection. 

periodic radio structure of PKS 0637-752. Regularly 
spaced trains of re-confinement shocks are routinely ob- 
served in hydrodynamic simulation s of AGN jets (e.g. 
lAlov et al.lll999l : lSaxton et al.ll2010( ). and synthetic radio 
maps based on numerical radiative transfer calculations 
indicate that simple luminous knots are associate d with 
the r e-confinement, or pinching shocks (S axton et al.l 
l2010f ). Analysis of the re-confinement process provides 
an estimate of the shock separation in terms of the jet 
power and external pressure. Let AR be the separation 
between re-confinement shocks, Qjct the jet power and 
Pcocoon the cocoon pressure. Analytic treatment of the 
re-confinement process, and numerical simulations of ax- 
isymmetric, hydrodynamic, relativistic jets indicate that, 
for a jet with Lorentz factor F = 5, 



AR w 65 



Qjc 



10^6 ergs/s 



1/2 



^cocoon 



10 



-11 



dyn/ c 



-1/2 



kpc , 

and is not se nsitive to the .let Lorent z factor or .iet 
opening ang le fKomissarov fc FalldHOOSl ) . iCroston et all 
|2005) find that the typical magnetic field strength in the 
lobes of high power FRII radio galaxies is in the order 
of 10 /i-G, and that the particle energy density is typ- 
ically a factor of a few greater than the magnetic en- 
ergy density. This suggests that the typical lobe pres- 
sure is in the order of 10~^^ dyn cm^^. Therefore, com- 
bining Equation [5] with the de-projected knot separa- 
tion (50 kpc < < 110 kpc) and an assumed value 
for the cocoon pressure (pcocoon ^ 10~^^ dyn cm"^), we 
find that the re-confinement shock interpretation implies 
Qjet« 10^6 ergs s-i. 

Hydrodynamic simulations by Sax ton et al.l (|2010f ) in- 
dicate that the ambient medium infiuences the distribu- 
tion of jet knots resulting from reconfinement shocks. 
In a constant pressure atmosphere, the jet knots re- 
main equally spaced along its length, but in an atmo- 
sphere with radially decreasing density, the knots become 
more closely spaced with distance along the jet. The 
nearly constant knot separation observed in PKS 0637- 
752 would therefore require the halo to have a core radius 
of hundreds of k pc, which is significa ntly larger than typ- 
ically observed ()Belsole et al.l |200"7[) . Furthermore, the 



re-confinement shock interpretation does not explain why 
the jet brightens in both the radio and X-ray bands at ap- 
proximately 8 arcseconds from the core. If the difference 
in jet brightness between the inner and outer jet were 
due to a difference in jet power, or a difference in cocoon 
pressure, then we may expect to see this reflected as a 
difference in knot separation. Instead, we see that the 
distance between knots remains approximately constant, 



independent of the jet brightness. 

3.1.2. Alternative Mechanisms 

A variety of mechanisms have be en proposed to explain 
regul a r patterns in AGN j e ts (e. g. iLapenta fc Kronberd 
l2005t IBalsara fc Normani [1991 IBahcah et al.l 119951 ). 
Most notably, Ke lvin-Helmholtz (KH) instabilities (e. g. 
iBirkinshawl (l990D have been successfully employed to 
describ e the emission patterns in pc-scale jets (see re- 
view bv lPeruc"holl2012f ). as w ell as the kpc-scale morphol- 
ogy of the jet in M 87 (e. g. iBicknell fc Begelman ,19961 
iLobanov et al.l[2003f ). A KH instability interpretation of 
the regular knot spacings in PKS 0637-752 would imply 
that the gross structure of the jet changes little over the 
region where this knot pattern is seen - in particular that 
the flow speed is constant - consistent with the IC / CMB 
model for the jet X-ray emission. We defer a more de- 
tailed investigation of the KH instability interpretation 
to an upcoming publication. 

3.2. Variation in the jet engine 

In the following sub-sections we consider models that 
invoke modulated activity of the jet engine to explain 
the semi-regular periodic appearance of the kpc-scale 
knots. To enable quantitative analysis, we must esti- 
mate the modulation timescale. Let us assume that the 
knots travel with Lorentz factor F, at an angle to the 
line of sight 6. The corresponding Doppler factor is 
S = [F(l — f3 cos9)]^^ . Accounting for projection and 
light travel time effects, the distance between the jet 
knots (D) is related to the apparent separation between 
the knots (Uapp) via 



D 



T5 sin 9 



(3) 



If the knot velocity /3 w 1 then V5 sin 9 k, /Jg^pp where 
/3app is the apparent transverse velocity of the kpc-scale 
knots. Hence, 



D 



app 



<D <D, 



app 



1 



Jcos0 



sin 9 



(4) 



The apparent transverse velocity of the pc-scale jet of 



PKS 0637-752 is 



'app 



13.3 ± 1.0 (|Edwards et al]|200l 

iLovell et al.ll2000( ) . We reasonably assume that the ap- 
parent transverse velocity of the kpc-knots is no greater 
than the apparent transverse velocity of the pc-scale 
knots, and therefore argue that /3app,max = 13.3. The 
one-sidedness of the kpc-scale jet implies that the ve- 
locity remains at least mildly relativistic on kpc-scales. 
Based on the current radio and X-ray images, the jet 
to counter-jet brightness ratio is > 60, implying that 
/3 > 0.5. As argued in Section [31 the jet viewing an- 
gle is 6* > 4°. These limits, combined with equation [H 
imply that if the knots are associated with moving jet 
components, 

0.6 kpc < £> < 60 kpc (5) 
and given 1 < /3 < 0.5, the inferred modulation period is 



2 X 10^ yr < r < 3 X 10^ yr. 



(6) 



The lower end of this range (r ~ 2 x 10'^ yrs) is applicable 
if the jet remains highly relativistic on kpc-scales, as 
suggested by the IC/CMB model of jet X-ray emission 
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(|Tavecchio et al.l 120000 . Given these constraints on the 
modulation timescale, we now consider two possible 
physical processes that could result in periodic modula- 
tion of the jet engine: accretion disk instabilities, and a 
binary supermassive black hole. 

3.2.1. Accretion disk instabilities 

In the context of accretion disk instabilities, we em- 
phasise a possible analogy between the jet of PKS 0637- 
752 and the quasi-periodic jet modula tion in the micro- 
quasa r GRS 1915-1-105 in the /3 state (|Fender fc Bellonil 
12004 and references therein). The variations in jet out- 
put in GRS 1915-1-105, with a characteristic timescale 
of ~ 30 minutes, are believed to result fr om limit cy- 
cle behaviour in an unstable accretion disk (IJaniuk et al.l 
[200l IJaniuk fc CzernvllMlt ITagger et al.ll2004Jl . Phvs- 
ical timescales are expected to scale linearly with the 
mass of the black hole, and therefore, the r ~ 30 minute 
oscillations in the jet of GRS 1915-M05 (A/bh « 14 Mq; 
IGreiner et al.| [2001) extrapolate to r ~ 2 x 10^ yrs fo r 
PKS 0637-752 (Mbh ~ 5 x 10^ Mq; ILiu et al.l[2006h : 
consistent with the central engine modulation timescale 
inferred from the periodic structure in the jet. 

It has previously been suggested that thermal-viscous 
instability of the accretion disk may play a role in shap- 
ing the morphology of arcsecond-scale AGN jets (e. g. 
ILin fc Shiel(SlT98alStaWarz et al.|[200l . Two modes of 
thermal-viscous instability potentially operate in accre- 
tion disks around massive compact objects: ionization 
instability and radiati on pressure instability (see e. g. 
IJaniuk fc CzCTml[20nL for a recent discussion). These 
instabilities are expected to result in quasi-periodic out- 
bursts due to modulation of the accretion rate, M . The 
ionization instability operates in the "partially ionized 
zone" in the outer disk, where hydrogen transitions from 
being neutral to ionized. The characteristic ionization 
instability timescale for sup ermassive black holes with 
Mrh ~ 10^ MfT, is > 10 ^ yrs (IMineshige fc Shieldsl [19901: 
IJaniuk fc Czernvl |2011[ ). which is longer than the out- 
burst timescale required to explain the periodic structure 
in the jet of PKS 0637-752. 

The radiation pressure instability operates in the in- 
ner disk where radiation pressure dominates over gas 
pressure, and results in short er timescale variability t han 
the ionization instability fsee lJaniuk fc Czernvll20lTI ). In 
the context of AGN, the radiation pressure instability 
has been proposed as an explanation for the intermit- 
tency of radio galaxies on short (~ 10^ yr) timescales, 
and in particular, the over-ab undance of young ra- 
dio sources in population studies (IRevnolds fc BegelmanI 
[19971 IMarecki et all 120031 ICzernv et al.ll2009D . At high 
accretion rates (near Eddington), and black hole masses 
Mbh ~ 10^ — 10^ M©, the characteristic timescale be- 
tween outbursts resulting fr om the radiation pressure in- 
stability is ~ lO'^ - lO'' yrs (jMerloni fc NavakshinI [20061: 
ICzernv et al.ll2009[ ) — consistent with the inferred mod- 
ulation period in PKS 0637-752. 

The interplay between the ionization and radia- 
tion pressure instability in AGN accretion disk s may 
affect the disk behavi our (jSiemiginowska et al.l 119961 : 
IJaniuk fc Czernvl 1201 ID . Variability over a range of 
timescales could potentially explain the variation in knot 
brightness between the inner and outer knots, as well as 



the variation in brightness between the knot and inter- 
knot regions. 

An alternative instability, the magnetic flooding 
accretion-ejection instability, has been proposed to ex- 
plain the /? state oscillations in GRS 1915+105, and may 
be applicable to AGN (T agger et al.ir2004V 

3.2.2. Binary supermassive black hole 

In this section we consider the possibility that the pe- 
riodic structure in the jet of PKS 0637-752 arises from 
variations in the accretion rate driven by a binary black 
hole. 

Supermassive black hole (SMBH) binaries are ex- 
pected to result directly from the merger of two massive 
gala xies, and may populate up to 10% of loca l galax- 
ies (jVolonteri et al. [ 2003; B urke-SpolaoH I2011D . How- 
ever, they are notoriously difhcult to detect; no defini- 
tive examples have been identified, although several po- 
tential binary precursor systems (i. e. dual active nu- 
clei in a single gal axy at large separation) have been 
confirmed (Komoss a et al.1 120031 : [Rodriguez et al.l [20061 : 
iFabbiano^t al.ii201lT) . Periodicities in morphology, mo- 
tion, and flux may arise from a bound supermassive 
black hole pair in a galaxy's centre (e .g. Ilguchi et al.l 
[20T0t [MacFadven fc Milosavlievic|[2008[ ). Excess emis- 
sion episodes correspond to a secondary SMBH's passage 
through disk-confined material (gas, stars) around the 
primary SMBH, inducing a heightened accretion rate at 
each traversal. Within this interpretation, the oscillatory 
signal indicates that the secondary SMBH's orbit is not 
in the same plane as the disk. Given a circular orbit, each 
observed knot may correspond to one disk passage event. 
However, if the orbit is elliptical, this model predicts that 
we are likely to observe double-peaked emission events 
corresponding to the initial and return traversal through 
the disk. The spacing of the sub-peaks will depend on 
the SMBH mass and ellipticity. OJ287, long theorised 
to be a SMBH binary, shows such structure in its light 
curve, emitting optical bursts at interval s of 11-12 years, 
with two sub-peaks per emission cycle ([Valtonen et al.l 
120081) . Future high-resolution observations of PKS 0637- 
752 may reveal additional structure in the jet, such as 
double peaked knots. 

Based on available information, we can estimate the 
properties of a putative binary system. The knot spac- 
ing corresponds either to P (for an elliptical orbit) or 
P/2 (for a circular orbit), where P is the orbital pe- 
riod of the binary system. This leads to periods of 
2 X 10'^ < P < 6 X 10^ years, and assuming a total sys- 
tem mass of ^ 5 X lO^M© (as estimated by ILiu et ahl 
[2006! based on luminosity and line-width), the corre- 
sponding orbital radius is 0.7 ^ a J; 30 pc. This range in 
orbital radii implies that the largest possible angular sep- 
aration of the binary lies in the range ~0.1-5 mas. How- 
ever, the current projected angular separation of the two 
SMBHs may be significantly smaller than the maximum 
angular separation, depending on the orbital orientation 
and phase. 

The binary black hole interpretation cannot explain 
why the jet brightens in both the radio and X-ray bands 
at approximately 8 arcseconds from the core, and an ad- 
ditional mechanism, such as variation in the jet power, 
must be invoked to explain the sudden increase in jet 
brightness at 8 arcseconds. 
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4. DISCUSSION AND CONCLUSIONS 

We have presented an 18 GHz ATCA map of the quasar 
PKS 0637-752, and identified a spectacular train of 9 
quasi-periodic knots extending 11 arcseconds along the 
jet. We sought to address the question "what physical 
process is responsible for the periodic knot structure?" , 
and considered two classes of model: (A) those that in- 
volve a static pattern through which the jet plasma trav- 
els, and (B) those that involve quasi-periodic modulation 
of the jet engine. One model that falls under class A is 
the re-confinement shock interpretation. If the knots are 
associated with re-confinement shocks, the observed knot 
separation implies that the jet kinetic power is approxi- 
mately 10"*^ ergs s^^. However, the constant knot separa- 
tion is not expected in a realistic external density profile. 
The re-confinement shock interpretation predicts a cor- 
relation between knot separation and jet kinetic power, 
which could be revealed in an imaging survey of a large 
sample of quasar jets showing regularly spaced knots. In 
this letter, we did not consider in detail alternative mech- 
anisms under class A such as Kelvin-Helmholtz instabil- 
ities. We defer a more detailed investigation of the KH 
instability interpretation to an upcoming publication. 

For models in Class B, the quasi-periodic structure in 
the jet is interpreted as arising from quasi-periodic mod- 
ulation of the central engine. The inferred modulation 
period is 2 x 10'^ yr < r < 3 x 10^ yr. The lower end of 
this range is applicable if the jet remains highly relativis- 



tic on kpc-scale s, as implied by the IC /CMB model of jet 
X-ray emission ()Tavecchio et al.ll2000D . This modulation 
timescale is consistent with the predicted radiation pres- 
sure instability timescale in the accretion disks of AGN. 
Indeed, we have drawn a direct comparison between the 
periodic structure in the jet of PKS 0637-752 and the 
quasi-periodic variation in jet activity seen in the mi- 
croquasar GRS 1915-1-105, believed to result from limit 
cycle behaviour in an unstable accretion disk. Finally, 
variations in the accretion rate may be driven by a sec- 
ondary black hole in orbit around a primary (i. e. more 
massive) black hole. We estimated the orbital radius of a 
putative binary system to be 0.7 ^ a < 30 pc, which cor- 
responds to a potentially-resolvable angular separation of 
^^0.1-5 mas. 
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